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Abstract 

We examine the effect of isospin-violating meson-nucleon couphng constants 
and of Ti-r] mixing on the binding-energy differences of mirror nuclei in a 
model that possesses no contribution from p-u: mixing. The ^He-^H binding- 
energy difference is computed in a nonrelativistic approach using a realistic 
wave function. We find the ^He-^H binding-energy difference very sensitive 
to the short-distance behavior of the nucleon-nucleon potential. We conclude 
that for the typically hard Bonn form factors such models can not account 
for the observed binding-energy difference in the three-nucleon system. For 
the medium-mass region (A=15-41) the binding-energy differences of mir- 
ror nuclei are computed using a relativistic mean-field approximation to the 
Walecka model. We obtain large binding-energy differences — of the order of 
several hundred keV — arising from the pseudoscalar sector. Two effects are 
primarily responsible for this new finding: a) the inclusion of isospin breaking 
in the pion-nucleon coupling constant, and b) the in- medium enhancement of 
the small components of the bound-state wave functions. We look for off-shell 
ambiguities in these results and find them to be large. 

PACS number(s): 21.10.Dr, 24.10.Jv, 13.75. Gx 
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I. INTRODUCTION 



Some of our most fundamental ideas about charge-symmetry-breaking (CSB) phenomena 
in nuclear physics [0,0 are currently being revisited. To a large extent this revision has 
been prompted by the suggestion of Goldman, Henderson, and Thomas of a suppressed 
contribution from p-uo mixing to the CSB component of the nucleon-nucleon (NN) interaction 
at small momentum transfers . Since then, many calculations, using a variety of theoretical 
approaches, have confirmed this suppression |^,|,^|^||,|[ . Although the issue continues to be 
controversial [0,^, the search for alternate CSB mechanisms has already started. Indeed, 
in a recent publication we have proposed isospin breaking in the meson-nucleon vertices as an 
additional source of isospin violation |jl2|. We have used isospin breaking in the pion-nucleon 
coupling constant to account for the large isospin violation in the pion-nucleon system 
reported recently by Gibbs, Ai, and Kaufmann [|1^]. Moreover, we have studied the impact of 
isospin violation in the vector-meson-nucleon coupling constants on the difference between 



the neutron and proton analyzing powers (AA) measured in elastic n-p scattering |T^. We 
have concluded that the magnitude of the resulting class IV CSB potential was consistent 



with that phenomenologically required by the experiment |T2 



However, other CSB observables, such as the scattering-length differences of the pp and 
nn systems and the binding-energy differences of mirror nuclei — the Okamoto-Nolen- 
Schiffer (ONS) anomaly |T7|JT8[] — are particularly insensitive to the class IV component of 



the CSB potential. Rather, these observables are driven by the class HI component. Indeed, 
several authors [^,^,^ have observed that the class HI component of the CSB potential — 
generated from "on-shell" p-uj mixing — dominates the binding-energy discrepancy of mirror 
nuclei and could explain most of the ONS anomaly 0. 

Motivated by the recent developments, we propose to study the ONS anomaly in a model 
which possesses no contribution from p-uj mixing; in our model, the resulting CSB potential 
is generated exclusively from tt-t] mixing and from isospin violation in the meson-nucleon 
coupling constants. To compute the '^He-^H binding-energy difference we employ a realistic 
three-body wave function obtained from the variational analysis of Ref. [22| using the Reid- 
soft core potential. To study the ONS anomaly in the A= 15-41 region we use a relativistic 
mean-field approximation to the Walecka model p3|]. Although the Walecka model has been 



used recently to compute binding-energy differences in mirror nuclei 1 , these studies 
have been limited to the vector sector: Coulomb corrections plus on-shell p-uj mixing. How- 
ever, since these contributions are dominated by the timelike part of the vector vertex (7°) 
they are not particularly sensitive to relativistic effects. In contrast, one expects important 
relativistic corrections in the pseudoscalar sector because of the structure of the (7^) ver- 
tex. Examining the role of relativistic effects and of isospin violation in the meson-nucleon 
coupling constants are two of the major components of the present investigation. 

Our paper has been organized as follows. In Sec. ^ we introduce the CSB potentials 
derived from meson mixing and from isospin violation in the meson-nucleon coupling con- 
stants. There, we review some of the ideas germane to our model. Sec. |T| is reserved to the 
analysis of the ONS anomaly in the three-nucleon system with special emphasis placed on the 
sensitivity of our results to the choice of hadronic form factors. Binding-energy differences 
of A=15-41 mirror nuclei are computed in Sec. ^ using various off-shell extrapolations of 
the NN potential. Finally, we present our conclusions in Sec. 0. 
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II. CHARGE-SYMMETRY-BREAKING POTENTIALS 



A. Vector-meson sector 



The CSB potentials which arise from isospin violation in the vector-meson-nucleon ver- 
tices and from p-uj mixing are given, respectively, by |1T|J^JT^ 



•^CSB 
•^CSB 
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Note that 9j^j<^^-i 9Njsfp, and /jvATp) are the isospin-conserving meson- nucleon coupling con- 
stants extracted from fits to two nucleon data and displayed in Table | using the Bonn B 
potential model Charge-symmetry breaking in the potentials is characterized by the 
appearance of three isospin-violating parameters: fi, /q, and {p\H\uj)] an isovector-tensor 
NNuj coupling constant, an isoscalar-tensor NNp coupling constant, and the p-u mixing 
amplitude, respectively. 

Our model is inspired by the notion of vector-meson dominance (VMD), so that the 
vector mesons couple to the appropriate isospin components of a conserved electromag- 
netic current. In models of this kind the p-u mixing amplitude necessarily vanishes at 
= . Moreover, our VMD assumption places important constraints on the form of 
the isospin-violating vector-meson-nucleon vertex at = 0. Specifically, isospin breaking 
can occur only in the tensor part of the vertex; the vector coupling is protected by gauge 
|l^,2^j2^. Thus, gauge invariance precludes the appearance of a class III — vector- 



mvariance 

vector — component in the CSB potential generated from vector-meson exchange. For the 
analysis of AA this fact is of no consequence, as class III CSB potentials do not contribute 
to np observables. Indeed, the class IV contribution from tu-meson exchange is able to fill 
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the role demanded by data; recall that this term is identical in structure and comparable in 
size to the one generated from on-shell p-uo mixing [12|. However, for observables in which 
the vector-vector component of V^gg was believed to be dominant — such as in the case of 
the ONS anomaly [p!9| , p0| , ^ — the absence of the corresponding term from V^gg and Vqo,^ 
could have important phenomenological consequences. This conclusion is unavoidable: it 
follows directly from gauge invariance; gauge invariance forces the isospin-violating vector 
coupling and the p-uj mixing amplitude to vanish at = 0. 



B. Pseudoscalar-meson sector 

The CSB potential which emerges from isospin violations in the pseudoscalar-meson 
sector is given by, 



(5) 
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(5) 



V^slB(QV(l)7'(2)k(l) + r.(2) 



'CSB 



(3a) 
(3b) 



(5) 

with V^gB made out of contributions arising from isospin violation in the pion-nucleon vertex 
and from tt-?] mixing |T^, respectively. 
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The isospin-conserving pion-nucleon coupling constant gj>^j^^ is given by the Bonn B 
potential value listed in Table |. For the NNr] coupling constant we have decided to 
use a value extracted from SU(3)-flavor symmetry rather than the Bonn-potential 
value of g'^j^^jAn = 2.25. Indeed, it is believed that the Bonn potential overestimates 
the coupling; a recent analysis based on ?7-photoproduction data suggests values as low 

as gNNr^/'^T^ < 



(see also Ref. 32 



For the r]N cutoff parameter we have simply 

assumed: 

violating pion-nucleon coupling constant {g^) and the n-rj mixing amphtude {{7i\H\ri)). CSB 
potentials derived from tt-t] mixing have been included in previous nonrelativistic analyses of 



the ONS anomaly |19,21,3C|. However, to our knowledge, CSB potentials arising from isospin 
violation in the meson-nucleon coupling constants have never been incorporated into these 
analyses, and only recently have they been used to study isospin violations in low-energy 
pion-nucleon scattering O . 



III. BINDING-ENERGY DIFFERENCE IN THE A=3 SYSTEM 

In this section we examine the effect of Tc-rj mixing and of isospin-violating meson-nucleon 
coupling constants on the ^He-^H binding-energy difference. Experimentally, this differ- 
ence is very precisely known and has the value of 764 keV. A nearly model-independent 
method for estimating the electromagnetic contribution to the binding-energy difference has 
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already been developed [^,0; ambiguities remain, however, on how to separate the meson- 
exchange- current component from the experimental charge form factors. These estimates 
suggest an electromagnetic contribution to the binding-energy of 693 ± 19 ± 5 keV, leaving 
71 ± 19 ± 5 keV to be explained by charge-symmetry-breaking effects [@,^. Note that a 
very similar result has been obtained recently by Ishikawa and Sasakawa using a different 



method |2T[. These authors have also computed the contribution to the binding-energy 
difference arising from (on-shell) p-uj mixing and from vr-r/ mixing, with the former play- 
ing the dominant role. Our model possesses no contribution from p-uj mixing. Hence, we 
would like to investigate the possible role of alternate CSB mechanisms in accounting for 
the non-electromagnetic contribution to the ^He-^H binding-energy difference: 

AE = {^He\VppfHe) - {^H\VnnfH) = -(^He\AV\^He) . (5) 

We shall compute the ^He-'^H binding-energy difference using the realistic three-body 
wave function obtained from the variational analysis of Numberg, Prosperi, and Pace p2| . 
These authors have written the nuclear Hamiltonian 
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2M 

in terms of the two intrinsic coordinates 



H=7^ipl + M)+J:V^^^ (6) 



(i<j) 





- ri - r^) , (7) 

and their corresponding conjugate momenta. Note that the indices 1 and 2 are reserved 
for the identical particles in the system; this is particularly convenient for the perturba- 
tive estimate of the ^He-^H binding-energy shift. The basic dynamical input is contained 
in the two-nucleon interaction {Vij) which was chosen to be the Reid soft-core potential 



3^ . Matrix elements of the Hamiltonian were computed using a truncated basis of 820 
harmonic-oscillator states and the ground-state wave function was obtained from direct ma- 
trix diagonalization. This wave function provides a satisfactory description of the charge 
form factors of ^He and ^H in the low-momentum-transfer region p2| , p5| . Moreover, a per- 
turbative estimate of the Coulomb shift in ^He yields 672 keV — a value consistent with that 



reported in Refs. p0| , |21 



For simplicity, the binding energy difference will be evaluated with the dominant ^Sq 
component of the trinucleon wave function. Note that for this case it is sufficient to evaluate 
the ^Sq component of the CSB potential: 

AV = VnnCSo) - VppCSo) . (8) 

The singlet CSB component of the NN potential is obtained by performing a nonrelativistic 
reduction of the various Lorentz structures given in Eqs. (|l|,|5) and by evaluating the meson 
propagators in the static limit ll|,0. We obtain (/f = f]\[]^p/g]^]\fp) 



qVM2 
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AV''(q) = -g 
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, qVM^ 
q2 _|_ 
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Note the appearance of the "relativistic-correction" factor q^/M^. This factor emerges from 
the allowed Lorentz structures present in the CSB potentials. Recall that — in contrast to 
previous analyses based on on-shell p-u mixing — no Lorentz vector structure [7'*(1)7^(2)] 
is allowed in our model. For the ir-r] mixing amplitude we have used the accepted value of 
{7r\H\ri) = —4200 MeV^ |^0|- For the isospin- violating meson- nucleon parameter we have 



adopted the nonrelativistic quark-model estimate of Refs. |T2 
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where m = 313 MeV is the average constituent quark mass and Am = 4.1 MeV is the 
down-up quark mass difference [0. Although these values are model dependent, they are 
insensitive to the spatial component of the nucleon wave function; they follow directly from 
its spin and flavor content. Moreover, the value of is compatible with other estimates 
available in the literature and, in particular, with the recent value reported by Henley 
and Meissner from QCD sum rules PE|] . 

To facilitate the interpretation of our results it is convenient to transform the above 
expressions into configuration space: 
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Note that in configuration space the ^He-^H binding-energy difference emerges from a simple 
overlap between the singlet component of the wave function [-R(a)] and the CSB potential 



AE = - a^da\R{a)\'^AV{V2a) = / I{a)da . 
Jo Jo 



(12) 



We will refer to these estimates as the point- coupling results, as hadronic form factor have yet 
to be introduced. We note, as a consequence of the relativistic factor q^/M^, the appearance 
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of a contact term [S{r)] in all single-meson exchanges (the contact term is absent from the n-r) 
mixing potential because of the additional meson propagator). For realistic wave functions — 
such as the one used here — the contact term does not contribute to the overlap as a result 
of the repulsive short-range correlations. Thus, we conclude that all single-meson-exchange 
potentials generate a repulsive AV, namely, a CSB potential that is attractive(repulsive) in 
the pp{nn) channel — in contrast to the experimental observation. Note that this conclusion 
depends exclusively on the sign of the isospin-violating meson-nucleon coupling constants; 
recall that in our model these constants are all positive as a result of the up quark being 
lighter than the down quark. 

This conclusion could change dramatically upon the inclusion of hadronic form factors. 
With form factors the contact term gets smeared over a distance of the order of A^^ (with 
A being the cutoff parameter) and could contribute to the overlap if this distance becomes 
larger than the size of the repulsive core. Most of our knowledge about hadronic form factor 
comes from phenomenological fits to a large body of two-nucleon data . These fits reveal 
hadronic form factors that are typically harder (A ^1.5 GeV) than the corresponding ones 
extracted from electron-nucleon scattering experiments (A~0.8 GeV). However, the subject 
of hadronic form factors is highly controversial. Indeed, (model-dependent) analyses of deep- 
inelastic-scattering (DIS) experiments suggest that the NNn form factor is considerably 
softer than AjVAr7r~l-5 GeV [3^,J^,4^. Moreover, Holinde and Thomas have shown recently 



that it is possible to obtain an adequate description of two-nucleon data with the softer 
pion-nucleon form factor (Aatat^ ^ 0.8 GeV) suggested by the DIS analyses These 



authors have also shown that the NNp form factor can also be made softer — but no softer 
than AAr^Tp — 1.2 GeV — without compromising the fit. Hence, we find it instructive to do the 
calculation of the ^He-^H binding-energy shift using, both, the hard form factors of Table | 
as well as the softer ones suggested by the analyses of Refs. [ PU| , ^P^^ . In the latter case 



we fix, for simplicity, both pseudoscalar cutoff parameters at ANNn{ri) = 0.8 GeV and both 
vector cutoff parameters at Annp{uj) = 1-3 GeV. We incorporate hadronic form factors into 
the calculation by modifying the point meson-nucleon coupling in the following way [p7| : 

9nN7t(v) ~^ 9NN7T(ri)(^^) = 9NN7T{r,)(^ + / ^1(t,)) ^ ^ (13a) 

dNNpiu.) ^?^7VpH(q') = ^?^7Vph(1 + qVAp(^))"' . (13b) 

In Table |T| we report the contribution to the ^He-^H binding-energy difference arising 
from n-r] mixing, vr— , p— , and cj— meson exchange, respectively. For comparison, we have 
also included the contribution arising from p-uj mixing, with the mixing amplitude fixed 
at its on-shell value: {p\H\uj) = —4520 MeV^ The first column of numbers displays 

the results in the point-coupling limit. As suggested, all single meson-exchange potentials 
give rise to a negative binding-energy shift; only n-i] mixing leads to a positive — albeit very 
small — contribution. The binding energy-shifts computed with the Bonn B hadronic form 
factors are reported in Table |T| under the heading of hard form factors. Evidently, the 
binding-energy shift is very sensitive to the short-distance behavior of the potential; the 
inclusion of form factors results in a 40% reduction relative to the point-coupling estimate. 
This reduction is best appreciated by plotting the integrand from which AE is obtained, 
that is I {a) in Eq. (|T^. In Fig. |I] we display the contribution from one-pion exchange to 
/(a); its contribution to AE for each estimate, which is the area under the curve, appears in 
parentheses next to its label. The solid line represents the point-coupling result. In this case 
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the contact term does not contribute and the integrand is negative-definite over the whole 
region. As the (hard) Bonn form factor is introduced (dashed fine) the smeared contact term 
dominates the short- distance region and cancels some of the negative contribution arising 
from the long-range term. It is then clear that for a sufficiently soft form factor the node 
in the potential — arising from the competition between the short- and long-range pieces — 
could be moved to a large enough distance as to make the integral positive; this is precisely 
what is observed when the softer form factor is introduced (dot-dashed line). The complete 
set of numbers with soft form factors has been collected on Table 0. For completeness we 
have also displayed in Fig. ^ the corresponding integrand for the omega meson. The same 
features are evident. In the present case, however, the integral is negative in all three cases 
because the vector form factor remains fairly hard. Nevertheless, we should mention that 
even for the case of (perhaps) unrealistically soft form factors, namely A7VAr7r(^) = 0.6 GeV 
and AArArp(^) = 0.85 GeV, our model can only account for +40 keV of the ^He-^H binding- 
energy difference; a result that is still one standard deviation away from the experimental 
value. 

In conclusion, we have shown that an analysis based on n-i] mixing and on isospin- 
violating meson-nucleon coupling constants — but one without p-u mixing — can not describe 
the '^He-^H binding-energy difference. However, we have established that this finding is 
extremely sensitive to the least understood feature of the NN potential, namely, its short- 
distance structure. Thus, it is difficult to draw any definite conclusion until an accurate 
determination of the cutoff parameters in the meson-nucleon form factors is made — a very 
difficult task indeed. 



IV. BINDING-ENERGY DIFFERENCES IN A=15 41 MIRROR NUCLEI 

In this section we examine the effect of meson mixing and of isospin violation in the 
meson-nucleon coupling constants on the binding-energy differences of medium-mass mirror 



nuclei using a relativistic mean-field approximation to the Walecka model [23|. The Walecka 
model is a renormalizable strong-coupling field theory of nucleons interacting via the ex- 
change of scalar and vector mesons. One of the great virtues of the model is its simplicity. 
Indeed, already at the mean-field level the Walecka model rivals some of the best available 
nonrelativistic calculations and provides a unified description of various nuclear properties, 
such as nuclear saturation, the spin-orbit force, and the density dependence of the nuclear 
interaction. The mean-field approximation is characterized by the existence of large Lorentz 
scalar and vector components that are responsible for a substantial in-medium enhancement 
of the small (lower) components of the single-particle wave functions; the so-called M* effect. 



A. Vector-meson sector 

In our model p-u mixing can not contribute to the ONS anomaly, at least at = 0. Yet, 
a calculation of binding-energy shifts using the on-shell value for the p-u mixing amplitude is 
still useful, as it provides a baseline against which alternate CSB mechanisms can be tested. 
Thus, we start this section by computing the impact of the dominant class III vector-vector 
component of Vqsb to the ONS anomaly. The purpose of this exercise is to establish the size 
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of the phenomenological gap that will have to be filled by the alternate CSB mechanisms. 
Note that we will not report the contribution to the binding-energy differences arising from 
the tensor-vector and tensor-tensor terms in Eq. ^ for spin-saturated nuclei their numerical 
impact is small (typically less than 10%). 

We compute the single-particle spectrum in a relativistic mean-field approximation to 



the Walecka model p3[ and incorporate CSB corrections in the Hartree-Fock approximation, 
i.e.. 



where 



+4E / (^^csB(q) pLt(q)p?/](q) , 
(iq 



-4E / 7ai^c^SB(q) 



(0)/ ■ ^ 



(14) 

(15a) 
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Note that the factor of four appearing in the above expressions emerges from the char- 



acteristic class III isospin factor 
orbitals in the Fermi sea, and we 
timelike-vector density 



Tz{l) + T^i^) . In Eq. ( [T3| ) the sum runs over occupied 
lave introduced the Fourier transform of the off-diagonal 



P(3a 



:q) = y'f/xW^(x)7V^'i-W«(x). 



(16) 



Several remarks are in order. First, we have ignored the contribution from the Dirac sea 
(negative-energy states) in the evaluation of the single-particle spectrum as well as in the 
CSB corrections to the energies. Second, for spherically symmetric nuclei only the timelike 
part [7°(1)7''(2)] of the potential contributes to the Hartree term; for static sources and 

Thus 



as 



fields, current conservation prohibits the appearance of a three-vector current p3 
far as the Hartree term is concerned, the CSB potential generated from p-u mixing behaves 
as a short-range Coulomb potential. Third, even though the spacelike-part [7(1) -7(2)] of 
the potential contributes to the Fock term, it is of little numerical significance and has been 
neglected. Finally, we suppress vector-meson production by ignoring retardation effects in 
the meson propagators. 

For spherically symmetric nuclei the calculation of the binding-energy shifts simplifies 
considerably. In this limit the eigenstates of the Dirac equation can be classified according to 
a generalized angular momentum k and can be written in a two component representation; 
i.e., 



1 



X 



fx) 



(^a = (a; m) = (n, k; m 



(17) 



The upper and lower components are expressed in terms of spin-spherical harmonics defined 
by 



3^«,m(x) = (x|/ijm) 
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if K < 0. 



(18) 
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The binding-energy shifts can now be easily evaluated. That is, 
where g = |q| and we have defined the following quantities 

occ 

^12(?) = E(2j. + i)pSi(?)pS(?), 

b 

occ 

= E(2J6 + 1) \{ja, -1/2; Jb, +l/2\L, 0)&iq) 



(19a) 
(19b) 

(20a) 
(20b) 



b\L 



Note that the Fourier transform of the off-diagonal timelike-vector density has become 



Pablil) 



dx 



ga{x)gb{x) + fa{x)fb{x) iiiqx) ; {la + lb + L = even) . 



(21) 



Using the Bonn B potential parameters of Table | |2^| we obtained the binding-energy 
differences reported in Table These values are in qualitative agreement with those 



reported by Blunden and Iqbal |19[ and, more recently, by Barreiro, Galeao, and Krein [25 



It is interesting to note that, in contrast to the long-range Coulomb potential, the Fock 
term generates a substantial exchange "correction" (of the order of 30-40%) to the direct 



Hartree contribution. Table [II also includes the remaining discrepancy between experiment 



and three theoretical calculations of the Coulomb displacement energy p6| , ^ . Note that 
the reported nonrelativistic discrepancies (Adme and Askii) differ — in some cases by a large 
fraction — from the recent relativistic estimates (Arel) of Koepf, Krein, and Barreiro ^6 . 



B. Pseudoscalar-meson sector: pseudoscalar NNtt coupling 

We shall now compute the binding-energy shifts that arise from the CSB potential given 
in Eq. (Rf) using a mean-field approximation to the Walecka model. Because of the assumed 



pseudoscalar vertex, parity precludes a Hartree correction to the energy [g^ . Thus, we only 
need to evaluate the Fock term. That is, 



(22) 



where we have introduced the Fourier transform of the off-diagonal pseudoscalar density 

p£^)(q) = I rfx Upi^h'e-^^-- W,(x) . (23) 

As in the previous section, the calculation simplifies considerably in the limit of spherically 
symmetric nuclei. In this limit, we obtain a binding-energy shift given by. 



(24) 
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where 



-^r^?) =E(2j6 + l)[0a,-l/2;j,,+l/2|L,0)pE^(g)] , (25) 

b;L 

and with 

pibL (9) = / dx\ga{x)fb{x) + fa{x)gh{x)\jL{qx) ] {la + h + L = odd) . (26) 

Jo L J 

Note that the pseudoscalar density is, as expected, hnear in the small components of the 
wave functions and, thus, highly sensitive to their in-medium enhancement. Moreover, 
the induced energy shifts are positive definite — as suggested by the ONS anomaly — for all 
single-particle states. 

The pseudoscalar contribution to the binding-energy difference of mirror nuclei is dis- 
played in Table |T|. We have listed separately the contributions arising from isospin violation 
in the pion-nucleon coupling constant and from ir-r] mixing — with the former accounting for 
about 60% of the total. Note that relativistic effects induce an additional enhancement in the 
binding-energy differences, relative to the nonrelativistic estimates of Blunden and Iqbal J19 



From Table |TTT| one observes that the binding-energy differences from the pseudoscalar sector 
are comparable — if not larger — than those computed originally with p-u mixing. Indeed, 
the pseudoscalar sector — alone — accounts for about 70-85% of the discrepancy between the 
nonrelativistic theory and experiment. 

To gauge the importance of relativity we have also computed the binding-energy differ- 
ences by assuming that the lower component of the wave function is determined from the 
free-space relation; the upper component remained unchanged, apart from a small normal- 



ization correction. These "nonrelativistic" values appear in parenthesis in Table |^ next to 
the uncorrected numbers. The relativistic corrections are, indeed, very significant. This is in 
sharp contrast to the relativistic enhancements — of less than 5% — observed for p-u mixing; 
recall that only the dominant vector- vector component has been included, which is relatively 
insensitive to the enhancement of the small components of the wave function. Perhaps the 
result that best captures the essence of the present analysis is the ^^Ca-^^K binding-energy 
difference (the state in Table [TV] ). Our baseline value for this difference is 65 keV; this 
represents our best attempt at reproducing the result by Blunden and Iqbal [|^, namely, a 
nonrelativistic calculation with = 0. However, introducing, both, relativistic corrections 
and isospin violation in the pion-nucleon coupling constant increases the baseline value to 
435 keV; this represents an enhancement of more than a factor of six. 



C. Pseudoscalar- meson sector: pseudovector NNir coupling 

The above results suggest a large contribution from the pseudoscalar sector to the ONS 
anomaly. Yet, a definite statement awaits, as these results could be sensitive to off-shell 
extrapolations. Off-shell ambiguities arise from the fact that different Lorentz structures for 
the A^A^-pseudoscalar-meson vertex — which are all equivalent on-shell — can generate vastly 
different predictions once the nucleon is off its mass shell. Indeed, a pseudoscalar (7^) vertex 
is very sensitive to the in-medium enhancement of the lower components of the nucleon wave 
function, while a pseudovector coupling [7^7^(p'— p)^/2M] is not. In a model in which chiral 
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symmetry is not explicitly realized — such as the one employed here — a pseudovector NNtt 
coupling seems to be preferred, as it is known to incorporate the correct low-energy pion- 
nucleon dynamics without sensitive cancellations MM. In contrast, a recent analysis of eta 



photoproduction data from the nucleon seems to favor a pseudoscalar NNr] coupling ||31|| . 

In this section we adopt a pseudovector representation for both {NNtc and NNt]) ver- 
tices. The advantage of such a choice is that it provides us with a lower bound for the 
binding-energy differences of mirror nuclei within the pseudoscalar sector. This, combined 
with the upper bound set by the pseudoscalar representation, should adequately reflect the 
sensitivity of our results to the various off-shell extrapolations. 

For the pseudovector case the angular- momentum algebra becomes slightly more com- 
plicated than before due to the explicit spin dependence of the operator. Recall that in the 
pseudoscalar case the spin dependence is implicit through the upper-to-lower coupling in- 
duced by the pseudoscalar vertex. The binding-energy shifts computed with a pseudovector 
NNtt coupling are now given by 



^) v?s'B(?)^r'(<!) 



(27) 



where 



b:J 



^f'^n^) = E(2J. + 1) 0a,-l/2;j,,+l/2|J,0)p£P(g) 



(28) 



The pseudovector density has been computed in the static [{p' — p)^ 0] limit and is given 
by 



„(PV) 
PabJ 



where 



r 

(g) = V \CjL{+Ka,+Kb) dxga{x)gbix)jL{qx) + 

L=J±1 ^ ^0 



CjLi-K.a,-K.b) dxfa{x)fb{x)iL{qx) ] {la + k + L = eveu) , (29) 

JO J 



(_l)(I^+J+l)/2 CL+ila- Ja)i2ja + 1) + {h " Jb)i2jb + 1) , if L = J + 1 



(2J + 1) [J- (la - Ja)(2ja + 1) " (h " Jb) (2j6 + 1) , if L = J - 1 . 

(30) 



Note that in the static limit a pseudovector vertex does not induce an off-diagonal coupling. 
Hence, relativistic effects, which appear as {f / gY corrections, are negligible. Moreover, the 
extra factor of that results from the pseudovector coupling is expressed in units of the 
free nucleon mass M — not in units of M* (see Eq. Thus, the binding-energy differences 
computed with a pseudovector vertex could be reproduced with a pseudoscalar vertex — only 
in the case in which the lower component of the wave function be determined from the free- 
space relation (see Table |V|). Thus, the square of the effective mass in the medium relative 
to its free space value (i.e., M*^/M^) should serve as an indicator of the sensitivity of the 
approach to the various off-shell extrapolations. Note that in a mean-field approximation 
to the Walecka model M*^/M^ is typically of the order of 40%. Thus, off-shell ambiguities 
are expected to be large (see Table 0). 
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The pseudovector contribution to the binding-energy difference of mirror nuclei is dis- 
played, alongside the pseudoscalar results, in Table 0. As suggested previously, these results 
are very close to those obtained in Table |V| using the free upper-to-lower ratio. More sig- 
nificantly, however, is the very dramatic reduction in the binding-energy differences relative 
to the pseudoscalar results. Indeed, while a calculation with a pseudoscalar vertex accounts 
for about 70-85% of the ONS anomaly, only 30% of it can be explained with a pseudovector 
vertex. It is important to stress that although some theoretical guidelines do exist, off-shell 
ambiguities might be difficult to resolve given that the contribution from the isospin- violating 
part of the potential to hadronic observables is small. 



V. CONCLUSIONS 

Binding-energy differences of mirror nuclei have been computed in a model that includes 
two sources of charge-symmetry-breaking in the NN potential; isoscalar-isovector mixing in 
the meson propagator and isospin violation in the meson-nucleon coupling constants. For 
the vector-meson sector we have used a VMD-inspired model. Thus, on the basis of this 
assumption — but little else — we concluded that neither p-uo mixing nor the dominant vector- 
vector component of the vector-meson-exchange potential can contribute to the binding- 
energy discrepancy between theory and experiment. In models of this sort this conclusion 
seems unavoidable — as gauge invariance forces, both, the p-uj mixing amplitude and the 
isospin-violating vector coupling to vanish at = 0. This is in contrast to previous analyses 
that employed the on-shell value of the p-oo mixing amplitude to explain the bulk of the ONS 
anomaly 0|T^j2^,2T]. Hence, in our model the anomaly must be explained exclusively in 



terms of Tc-r] mixing and of isospin violation in the pion-, rho-, and omega- nucleon coupling 
constants. 

To compute the ^He-^H binding-energy difference a nonrelativistic estimate of the ^5*0 
component of the CSB potential was obtained. In our model all sources of CSB generate 
a ^5*0 interaction with a long-range part that is attractive in the pp channel and repulsive 
in the nn channel. This result is robust, as it only depends on the sign of the tt-t] mixing 
amplitude and on the sign of the isospin-violating meson-nucleon coupling constants; recall 
that in our model the latter is determined from the up-down quark mass difference, which 
we assume to be negative. The ^He-^H binding-energy difference, however, was seen to be 
very sensitive to the short- distance behavior of the potential. Indeed, the binding-energy 
difference emerged from a delicate interplay between the short- and long-range parts of the 
potential. This interplay introduces a node in the CSB potential which can be shifted to 
a large enough distance as to make the binding-energy difference almost consistent with 
experiment. We stress that the (almost) agreement with experiment could only be achieved 
at the expense of introducing unrealistically soft form factors. 

For the medium-mass (A=15-41) region we computed the binding-energy difference of 
mirror nuclei in a relativistic mean-field approximation to the Walecka model. Due to the 
larger size of the system, relative to ^He, our results were fairly insensitive to the short- 
distance structure of the potential. Moreover, explicit spin-dependent effects are known to 
be small for these spin-saturated nuclei. As a consequence, the numerical impact from the 
vector-meson sector on the ONS anomaly was seen to be small. We proposed pseudoscalar- 
meson exchange as an alternate mechanism to explain the ONS anomaly in the medium-mass 
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region. We computed binding-energy differences that were comparable — and in most cases 
larger — than those obtained from previous estimates using on-shell p-uo mixing. We have 
concluded that the pseudoscalar sector — alone — could explain about 70-85% of the ONS 
anomaly. Two effects, ignored until now, were responsible for these findings: a) the inclusion 
of isospin violation in the pion-nucleon coupling constant, and b) relativity, which enhances 
the small components of the bound-state wave functions relative to their free-space value. 
We reached these conclusions by assuming a pseudoscalar coupling for, both, the NNn and 
NNt] vertices. Yet other choices — all of them equivalent on-shell — are possible; such as 
pseudovector coupling. There is no definite prescription on how to take the CSB potential 
off-shell. Moreover, this issue is complicated further by the fact that chiral symmetry favors a 
pseudovector representation for the NNtt vertex while a recent analysis of ?7-photoproduction 
data seems to suggest a pseudoscalar NNrj coupling. Hence, to gauge the sensitivity of our 
results to the various off-shell extrapolations we repeated the calculation by adopting a 
pseudovector representation for both vertices. We concluded that the off-shell ambiguities 
were large; while 70-85% of the anomaly could be explained with a pseudoscalar coupling 
only 30% of it could be accounted for with a pseudovector vertex. 

In conclusion, we have used a VMD-inspired model to examine the effect of isospin- 
violating meson-nucleon coupling constants and of tt-t] mixing on the binding-energy differ- 
ence of mirror nuclei. We could account for most of the ^He-^H binding-energy difference — 
provided very soft form factors were assumed. Moreover, we could explain the ONS anomaly 
in the medium-mass region — provided a pseudoscalar representation was adopted. Based on 
these findings we expect that the Okamoto-Nolen-Schiffer anomaly will remain — even after 
more than a quarter of a century — the source of considerable theoretical debate. 
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FIGURES 



FIG. 1. The integrand for the ^He-^H binding-energy difference arising from one-pion exchange 
for three different choices of form factors: Atvattt oo (soUd Une), A^nnw = 1-7 GcV (dashed hne) 
Aattvtt = 0.8 GeV (dot-dashed hne). The one-pion exchange contribution to AE, which is the area 
under the appropriate curve, appears in parentheses next to its label. 

FIG. 2. The integrand for the ^He-^H binding-energy difference arising from one-omega ex- 
change for three different choices of form factors: Aatatoi oo (solid line), Annuj = 1-85 GeV 
(dashed line) Ajvato; = 1-3 GeV (dot-dashed line). The one-omega exchange contribution to AE, 
which is the area under the appropriate curve, appears in parentheses next to its label. 
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TABLES 



TABLE L Meson masses, coupling constants, tensor-to-vector ratio, and cutoff parameters of 



the Bonn B potential. 

Meson Mass (MeV) c/V4vr f/g A (MeV) 

~ 138 14.21 — 1700 

r? 549 0.90 — 1700 

p 769 0.42 6.1 1850 

u 783 11.13 0.0 1850 



TABLE IL 


Contribution to the ^He-'^ 


H binding-energy difference (in 


keV) arising from tt-t? 


mixing, vr— , p- 


-, and w— meson exchange 


for three different values of the cutoff parameters; see 


text for details. The experimental value is 


AE' = 71 ± 19 lb 5 keV. For comparison, we also include 


the contribution from on-shell p-u) mixing 








Point coupling 


Hard form factors 


Soft form factors 


TT-r] 


+7.4 


+ 13.6 


+22.2 


TT 


-27.9 


-16.2 


+13.0 


P 


-49.4 


-37.3 


-20.4 


LO 


-53.8 


-40.9 


-22.4 


Total 


-123.7 


-80.8 


-7.6 



p-u +57.9 +77.9 +87.3 



TABLE III. Contribution to the binding-energy differences of mirror nuclei (in keV) arising 
from isospin violation in the pion-nucleon coupling constant, vr-r/ mixing, their sum, and on-shell 
p-uj mixing. Also included are the remaining differences (A) between experiment and the Coulomb 
energy computed in three different models [p^, 43 1. 





Pseudoscalar 




Rho-ome^ 


;a 


ONS Anomaly 


A 


State 


^Qnn^ vr-r/ 


Total 


Hartree 


Fock 


Total 


Adme 


AsKII Arel 


15 




201 


145 


346 


379 


-149 


230 


380 


290 


95 


17 


1^5/2 


83 


56 


139 


246 


-91 


155 


300 


190 


79 


27 


^"5/2 


223 


148 


371 


445 


-176 


269 


480 


490 


551 


29 


2Sl/2 


165 


108 


273 


305 


-91 


214 


290 


240 


34 


31 


^■^1/2 


202 


134 


336 


355 


-119 


236 


540 


560 


39 


33 


14/2 


159 


113 


272 


363 


-121 


242 


360 


280 


213 


39 


^"3/2 


255 


180 


435 


455 


-169 


286 


540 


430 


404 


41 


1/7/2 


133 


89 


222 


355 


-129 


226 


440 


350 


404 
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TABLE IV. Contribution to the binding-energy differences of mirror nuclei (in keV) arising 
from isospin violation in the pion-nucleon coupling constant, ir-rj mixing, their sum, and on-shell p-LO 
mixing. Included in parenthesis are the corresponding quantities computed with lower components 
generated from the free-space relation. 



A 


State 






Total 


Rho-omega 


39 


1^5/2 


205 (103) 


142 (73) 


347 (176) 


294 (291) 




^*l/2 


233 (84) 


158 (58) 


391 (142) 


283 (271) 






255 (92) 


180 (65) 


435 (157) 


286 (275) 


41 


1/7/2 


133 (78) 


89 (53) 


222 (131) 


226 (224) 



TABLE V. Contribution to the binding-energy differences of mirror nuclei (in keV) arising 
from isospin violation in the pion-nucleon coupling constant, Tr-ry mixing, and their sum. The 
first(second) set of numbers were computed using a pseudoscalar(pseudovector) NNir vertex. 





Pseudoscalar 


Pseudovector 


A 


State 




TT-r] 


Total 




TT-rj 


Total 


15 




201 


145 


346 


67 


49 


116 


17 


lc^5/2 


83 


56 


139 


52 


36 


88 


27 


1472 


223 


148 


371 


109 


74 


183 


29 


2si/2 


165 


108 


273 


54 


36 


90 


31 


2s-^ 

^*l/2 


202 


134 


336 


62 


43 


105 


33 


14/2 


159 


113 


272 


49 


35 


84 


39 


-^"3/2 


255 


180 


435 


83 


59 


142 


41 


1/7/2 


133 


89 


222 


74 


51 


125 
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